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ABSTRACT 

We consider the general problem of the tidal capture or circularisation from 
large eccentricity of a uniformly rotating object. We extend the self-adjoint 
formalism introduced in our recent paper(Papaloizou & Ivanov 2005 hereafter 
PI)) to derive general expressions for the energy and angular momentum trans- 
fered when the planet or a star passes through periastron in a parabolic or 
highly eccentric orbit around a central mass. These can be used without mak¬ 
ing a low frequency approximation as was done in PI. We show how these 
can be adapted to the low frequency limit in which only inertial modes con¬ 
tribute for baratropic planet models. In order to make quantitative estimates, 
we calculate the inertial mode eigenspectrum for planet models of one and five 
Jupiter masses Mj, without a solid core, with different radii corresponding 
to different ages. The spectra are found in general to be more complex than 
of a polytrope with index n = 1.5, considered in PI, because of the existence 
global modes associated with the transition from molecular to metallic Hy¬ 
drogen. Nonetheless the main tidal response is still found to be determined 
by two global modes which have polytropic counterparts. 

These also determine the uniform angular velocity in a state of pseudo 
synchronisation, for which the angular momentum transfered during an en¬ 
counter is zero. This is found to be close to 1.55 times the circular orbit 
angular velocity at periastron for all models considered. This is in contrast to 
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the situation when only the / mode is considered (Ivanov & Papaloizou 2004, 
hereafter IP) and the equilibrium angular velocity is found to be much larger. 

We consider the multi-passage problem when there is no dissipation find¬ 
ing that stochastic instability resulting in the stochastic gain of inertial mode 
energy over many periastron passages occurs under similar conditions to those 
already found by IP for the / modes. We find that this requires circularisation 
to start with a semi-major axis exceeding ^ 30AC7, for final periods of ^ 3 
days, reducing to ~ I — 2AU for final periods ~ 1.2 days. 

Finally we apply our calculations of the energy transfer during a perias¬ 
tron passage to the problem of the tidal circularisation of the orbits of the 
extrasolar planets in a state of pseudo synchronisation, expected because of 
the relatively small inertia of the planet, and find that inertial mode excitation 
dominates the tidal interaction for IMj planets that start with semi- major 
axes less than lOAU and end up on circular orbits with final period in the 
4 — 6 day range. It is potentially able to account for initial circularisation up 
to a final 6 day period within a few Gyr But in the case of 5Mj oscillation 
modes excited in the star are more important. 

Key words: hydrodynamics; stars: oscillations, binaries, rotation; planetary 
systems: formation 


1 INTRODUCTION 

The process of tidal capture, whereby a body in initially unbound orbit has a close approach 
to another gravitating mass, and through the excitation of internal modes of oscillation, 
becomes subsequently bound is of general importance in astrophysics. It is believed to play 
a role in binary formation in globular clusters (eg. Press & Teukolsky 1977, hereafter PT, 
and references therein) and to play a role in the tidal interactions of stars in galactic centres. 
A related problem is the circularisation of an orbit from large eccentricity which may occur 
after a tidal capture. In this situation, each periastron passage is very similar to a close 
encounter in a weakly unbound orbit that causes a change in the energy contained within 
internal modes of oscillation at the expense of that in the orbit. 

In order to estimate capture probabilities and circularisation time scales the excitation of 
normal modes of oscillation by a perturbing potential must be considered. In general modes 
associated with spherical stars, or small perturbations of them, have been considered (eg. PT, 
IP and references therein). However, if the tidally perturbed body has multiple encounters. 


* E-mail:P.Ivanov@damtp.cam.ac.uk (PBI) J.C.B.Papaloizou@damtp.cam.ac.uk (JCBP) 
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Dynamic tides in rotating objects 3 
angular momentum transfer should lead it to rotate at an angular velocity related to that at 
periastron. In this case modes of oscillation with periods comparable to the rotation period 
need to be considered (eg. Papaloizou & Pringle 1978). 

An approach for dealing with the problem of calculating tidal phenomena while taking 
into account the part of the oscillation spectrum with eigen frequencies comparable to the 
rotation frequency was recently proposed by PI. They adopted a self-adjoint formalism for 
calculating the tidal response that used the density perturbation to describe the motions 
within the star and they made a low frequency approximation that the rotation and mode 
frequencies were very much less than the inverse sound crossing time. They also evaluated 
the inertial mode spectrum for a polytrope with index n = 3/2 and showed that inertial 
modes dominated the response at low frequencies or large periastron distances. 

Given the importance of rotation, the purpose of this paper is to develop further the 
work begun in PI and to use the results to obtain time scales for tidal phenomena such as 
circularisation of the orbit starting from high eccentricity. 

We begin by noting that the self-adjoint formalism can be extended to general linear 
displacements without making a low frequency approximation. In this form, which uses the 
Lagrangian displacement to describe the perturbations it has wider applicability than the 
formalism presented in PI which made a low frequency approximation. We also wish to 
evaluate the expressions obtained for the energy and angular momentum exchanged during 
close periastron passages for realistic planet models. Accordingly we calculate sequences of 
models of one Jupiter mass Mj and 5Mj. Each sequence has a range of radii in order to take 
account of the fact that the planet radius varies with age. For simplicity we do not include a 
solid core noting that according to Saumon & Guillot (2004) observations of Jupiter are not 
inconsistent with this. We calculate the oscillation spectra together with required overlap 
integrals for the models using a set of basis or trial functions. 

We then use the determined eigenmode properties for these models to estimate tidal 
phenomena such as the attainment of pseudo synchronism for which the angular momentum 
transfer at periastron passage is zero. Such a state is expected to be attained, on account 
of the small inertia of the planet compared to that of the orbit, with subsequent energy ex¬ 
changes occurring with the planet rotating at an equilibrium angular velocity. The inclusion 
of inertial modes is important for determining this angular velocity which turns out to be 
close to that of the orbit at periastron. Galculations in IP incorporating only the / modes 
gave a much higher value. 
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Having calculated the equilibrium angular velocity we are able to consider multiple en¬ 
counters and determine conditions for the onset of stochastic instability, that results in the 
stochastic gain of inertial mode energy over many periastron passages when there is no dissi¬ 
pation (Kochanek 1992, Mardling 1994a,b). Then we make estimates of circularisation time 
scales for the two different planet masses taking into account the evolutionary variation of 
the planet radius with age. 

The plan of the paper is as follows. In section 2 we formulate the new self-adjoint for¬ 
malism for the problem of linear perturbations proposed in our recent paper(PI) using the 
Lagrangian displacement to describe the motions and then use it to derive general expres¬ 
sions for the energy and angular momentum transfered when a rotating planet or a star 
passes through periastron in a parabolic or highly eccentric orbit around a central mass 
without making a low frequency approximation as in PI. It can thus be applied to general 
situations in which p, /, g and inertial modes are signihcant. In section 3 we show that the 
general formalism can be reduced to that of PI in an appropriate low frequency limit. In this 
case only the inertial modes contribute to the tidal response for the baratropic stellar mod¬ 
els we consider. In section 4 we discuss the tidal energy and angular momentum transfered 
through inertial modes giving very simple expressions that apply when the planet or star 
rotates in a state of pseudo synchronisation for which the angular momentum transfered in 
an encounter is essentially zero. We apply these to the multi-passage problem when there 
is no dissipation hnding conditions for the occurrence of the stochastic instability. We hnd 
that this requires the circularisation process to start with a rather large semi-major axis 
~ 30AU, for hnal periods of ~ 3 days. However, this could be reduced to ~ 1 — 2AU for the 
shortest observed hnal periods ~ 1.2 days. 

In section 5 we calculate the eigenspectrum of the inertial modes for the planet models 
which have been obtained using a realistic equation of state. We hnd that the counterparts 
of the two global modes found in a polytrope with n = 1.5 exist in these cases but that 
there are additional global modes associated with the transition from molecular to metallic 
Hydrogen. We go on to apply the results to the problem of tidal circularisation of the orbits 
of the extrasolar planets in sections 6 and 7 hnding that inertial modes dominate the tidal 
interaction for Mj planets ending up on circular orbits with hnal period in the 4 — 6 day 
range and are potentially able to account for the initial circularisation at a 6 day period 
within a few Gyr in that case. However, in the case of 5Mj oscillation modes excited in the 
star become more important. Finally in section 8 we discuss our results. 
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2 FORMAL SOLUTION OF THE TIDAL PROBLEM 

In this Section we further develop the self adjoint formalism describing the pulsations of 
a uniformly rotating fully convective object (a planet or a star), referred to hereafter as a 
planet, introduced in our previous Paper (Papaloizou & Ivanov 2005, PI). We give general 
expressions for the energy and angular momentum transfer from the orbit to the planet as 
a result of a close encounter or fly by around a gravitating centre. Here, contrary to PI, we 
do not make the assumption that only low frequency modes respond effectively to the tidal 
forcing and the results obtained are valid for the excitation of any pulsational mode of a 
rotating planet. 

The fundamental quantity used in our analysis is the Lagrangian displacement vector 
By dehnition, this is real. We adopt a cylindrical coordinate system {w, 0, z) with origin at 
the centre of mass of the planet. 

We make use of the Fourier transforms of the displacement vector and any other pertur¬ 
bation quantity, say Q, over the azimuthal angle 0 and time t in the form 

/ -I-OO _ \ 

dcrQm exp(—iaf) + cc j , ( 1 ) 

-OO / 

where the sum is over m = 0 and 2 and cc denotes the complex conjugate. The reality of Q 
implies that the Fourier transform, indicated by tilde satishes Qm(o') = Q-rni~^)- 

The inner products of two complex scalars Yi, Y 2 , and vectors and r ]2 are dehned as 

(Yi|y2) = J zudzudzp{Y^Y 2 ) {'ni\'n2) = J wdwdzp{rjl ■ 772), ( 2 ) 

where (77* ■ 77 ) is the scalar product, p is the density, and * stands for the complex conjugate. 

2.1 Perturbed equations of motion 

We assume that, viewed from an inertial frame, the planet is rotating with uniform angular 
velocity H directed perpendicular to the orbital plane. In this situation the hydrodynamic 
equations for the perturbed quantities take the simplest form in the rotating frame which 
we accordingly adopt. 

Since the planet is fully convective, the entropy per unit of mass of the planetary gas 
remains approximately the same over the volume of the planet and the pressure P can 
be considered as a function of density p only, P = P{p). A perturbation of the planet may 

^ Note that this condition can be broken at the radius corresponding to phase transition of hydrogen from metallic to molecular 
form, see below. 
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be considered as adiabatic, and therefore the same fnnctional dependence P = P{p) holds 
dnring pertnrbation as well. This condition is often referred to as the baratropic condition. 
In the baratropic approximation the linearised Navier-Stokes eqnations take the form (see 
PI) 

(3) 




where 

W = cIp + + (4) 

p! is the density pertnrbation, Cg is the adiabatic sonnd speed, is the stellar gravitational 
potential arising from the perturbations and is the external forcing tidal potential. The 
viscous force per unit of mass is f^. In Cartesian coordinates it can be written in component 
form as 


2 • 


/; = -(i/rtc/ = (i/p){pKCf ,+A - 


(5) 


where the comma denotes differentiation and u is the kinematic viscosity. The Greek indices 
enumerate spatial coordinate directions and the usual summation convention over these is 
adopted. 

Note that the centrifugal term is absent in equation being formally incorporated 
into the potential governing the static equilibrium of the unperturbed star. The convective 
derivative ^ ^ as there is no unperturbed motion in the rotating frame. 

The linearised continuity equation is 

p' = -V-(pO, (6) 

and the linearised Poisson equation is 

AT*”* = AnGp. (7) 


Provided that expressions for the density, sound speed and kinematic viscosity are spec- 
ihed in some unperturbed model of the planet, the set of equations (EHID is complete. 
Without the viscous term on the right hand side of equation © these equations coincide 
with the analogous equations used by PI. 

The perturbed potential T*”* can be formally expressed in terms of the density pertur¬ 
bation with help of the inverse of the Laplace operator, A“*, and the density perturbation 
is expressed in terms of the displacement vector with help of the continuity equation. We 
have 
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=-AttGA-^V ■ {p$,). (8) 

Now we express W in terms of ^ with help of equations ©, ® and ®. We substitute 
the result into equation thus obtaining an equation containing only the displacement 
vector ^ as a dynamical variable. After making Fourier transform of this equation we get 

“ f = ^m, (9) 

where the operators B and C are defined through their action on the vector ^ as 

Bi = -2in X (10) 

and 

C^ =-V((^ + 47rGA-^)V-(p^)) (11) 

P 

It can be easily checked that the so defined operators B and C are self adjoint in the sense of 
the normal product Q. Also for the models we consider with weak self-gravity, the operator 
C may be taken to be positive dehnite for displacements with a non zero density perturbation 
somewhere and non negative for arbitrary displacements including those with zero density 
perturbation everywhere. 

The source term Sm is determined by the external forcing potential 

5 ™ = ( 12 ) 

2.2 Self adjoint form of the perturbed equation of motion and its formal 
solution 

The presence of the viscous force leads to dissipation of energy of the perturbed motion 
on a dissipation time scale. For global disturbances in astrophysical systems this time scale 
is typically very large compared to a time scale characterising the tidal forcing. In this 
situation we can treat the action of the viscous force as a perturbation. 

Neglecting the contribution of the viscous term, we see equation ® has an analogous 
form to one discussed by PI for low frequency perturbations. There it was shown that this 
form led to a standard type of eigen value problem associated with a new self adjoint operator 
7i. In order to bring equation © to this standard type, let us formally consider the square 
root of the operator C, defined by the condition: C = This can be done on 

account of the non negativity of C. Further, let us introduce a six dimensional vector Z with 
components 
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Zi = <^L. (13) 

It follows from equation m that the eigenvalue problem can be written in the standard form 
aZ = nZ + S, (14) 


where the operator 7i has a matrix structure 

( B \ 

[ 0 ) 

and the source six dimensional vector S has the components 
Si = Sm, S 2 = 0. (16) 



It follows from equation ()15|1 that the operator T-C has the required self adjoint form. Now 
we can look for solution to equation m in the form: 

Z = Y,akZ\ (17) 

k 

where Z^ are the eigenfunctions of Ti which satisfy 

akZ^ = nZ^. (18) 


The eigen frequencies ak are necessarily real. Equation m is equivalent to equation m 
with the source term Sm = 0. 

Note that the eigen spectrum of or a part of it can be continuous. It this case the 
summation in (HHD must be replaced by integration with an appropriate measure. 

The eigen functions Z^ are orthogonal in the sense of the inner product 

< Z’‘\Z‘ >= (ZflZ'i) + (Z‘|Z') = (19) 


where = Zj/trr 

Substituting (HID into equation dn and using (USD we obtain 

^ ^ < Zk\S > ^ alSk ^ 20 ) 

^ Nk{a + ia^ - ak) Nk{a + ia^ - ak) ’ 

where we add a small imaginary part to the frequency a with Uj, > 0 in potentially resonant 
denominators to account for dissipation. As we show below, for modes for which the effect 
of viscosity can be treated as a small perturbation, an explicit expression for cTj, can be 
obtained by taking into account the dissipation determined by the viscous force f^. 

N, =< Z,\Z, >= a|(?^|?^) + (21) 


is the norm. 
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The coefficients Sk determine decomposition of the sonrce vector over the eigen fnnctions, 

St, (22) 


- _„<Z^ _^4St 

,n 2^ N. 1 


k ^ k 

where snmmation is performed over all modes with a particular value of m. They can be 
written in the form 

S, = {^,\~S^) = {pX'). (23) 

where is the density perturbation associated with the eigenvector = —Vp^^, and 

to obtain the last equality, we integrate by parts. 

It is important to note that the coefficients Sk, the eigen frequencies ak and the associated 
displacement vector obey a certain identity relation 

(24) 


0 


t Nt 

which follows from equations (USD, (|22D and the orthogonality relation (USD for the eigen 
vectors. This relation is also required to ensure that remains bounded as a —0 (see 
below). 

The solution for the displacement vector follows from equations (USD, (ED and (ED 
-^ISk ^ V- (^kSk 


e =y — 

^ Nka{a + zo-j, - ak) 


Nk{a+ ia^, - ak) 


^ki 


(25) 


where we use the relation ()24|1 to obtain the last equality. 

The displacement vector ^ is obtained from its Fourier’s transform by integration 
over a and summation over m 


^ = M- 


^kSk 




(26) 


^ Nk{a + iay - ak) 

This expression can be simplihed in the limit f —> cxo at which only the poles in the expression 
in the braces contribute signihcantly to the integral over a. In this limit we may, therefore, 
write 


« = 2" + C.C.). 

m,k 


(27) 


2.3 Energy and angnlar momentnm transfer 

It can be easily shown from equation dSD or equation m that the expression for the canonical 


energy of the perturbations is given by 

JV Z rm 


(28) 
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where we integrate over the volume of the planet. 

We evaluate the expression (12HD after a close encounter, t > 0. Assuming that the time 
after closest approach is much larger than the characteristic time scale of the encounter, but 
nevertheless much smaller than the characteristic viscous time scale ~ , we can use 

equation (EZD setting (Ty = 0. Substituting the expression for ^ and its time derivative into 
and taking into account equation m we get 
<ylSuSl 


E, = 87r^Y. 

m,k 


Nt 


(29) 


where we use the orthogonality relation (P|) and the expression for the norm (HI. 

The rate of viscous dissipation of energy is given by the standard expression 

(30) 


J + O - 5(C)d- 


An explicit expression for the amount of energy dissipated by viscosity after the fly by, 
AEi, = dtE^, can be obtained from dSHD substituting into dsni) and integrating over 
time the resulting expression. We have 

<^k(^lSk 


AE^ = Stt^ / daa^ Y ——-- YkSiaki, 

' ^^lNkNi{(T - (Jk+l(Ty){(T - (Ti-KT^) 


(31) 


where we use the well known relation / dte''^^~^ > = 2'n8{a — cr), and 

a,, = I wd,^dz{l;lMQjr + (CT) - 5e,»(?w)*}' (32) 

As we have mentioned above, we assume that viscosity is very small, and the imaginary part 
of the frequency is much smaller than the corresponding real part. In this case, one can see 
from equation m that the leading contribution to AEy is determined by diagonal terms 
in the summation series with k = 1. Also, only a part of the integral over a very close to the 
resonance, a ~ <7^, gives a signihcant contribution. We obtain 




da 


N? 

m,k '^^k 


^k.k 


Svr^ _ aiSkSt 


E 


m,k ^k 


~^k,k • 


(33) 


i<z-aty‘ + al) 

The amount of energy transfered by tides into the modes must be equal to the amount of 
energy dissipated by viscosity: Ec = AEy. From this condition and equation (12ni) and dnsi) 


we easily obtain an explicit expression for the inverse dissipation time 


(Ty = 


2 

^k^k,k 

Nt 


(34) 


Similar arguments can be used for calculations of the amount of the angular momentum 
transferred to the modes after the fly by. However, the expression for the angular momentum 
transfer directly follows from the general relation between canonical energy and angular 
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momentum corresponding to a particular mode with given values of k and m (eg. Friedman 
and Schutz 1977) 

Lc = mEc/ak- (35) 

Accordingly, we have from equation {a 

(36) 

2 k 


The expression for the energy of the modes in the inertial frame, 77/, follows directly 
from equations and (jSi 



Note that the modal energy in the inertial frame can be negative. 


3 RESPONSE OF LOW FREQUENCY MODES TO THE TIDAL FORCING 

The formalism developed in the previous Section can be applied for any pulsational mode of a 
rotating planet. It follows from equations (El, El and El that in order to hnd the transfer 
of energy and angular momentum from the orbit to the modes we should be able to hnd 
solutions (eigen functions and eigen frequencies ) of equation (fTH|l describing free oscillations 
of the planet. This can be done numerically. However, in general, it is rather difficult to use 
numerical methods directly in solving equation El without any further simplihcations. It 
turns out that these simplihcations can be naturally introduced into the problem on hand. 
Assuming that the periastron distance Vp is sufficiently large, there is a small parameter in the 
problem, namely the ratio of the typical frequency associated with the periastron passage, 
~ typical frequency associated with the planet, We 

would expect that the angular momentum transfer between the orbit and the modes leads to 
a so called state of pseudo synchronisation where H ~ flp, and therefore the angular velocity 
of the planet can be, typically, much smaller than the characteristic frequency H*. In such a 
situation all modes of pulsation present in a convective planet can be divided in two distinct 
groups: high frequency modes with eigen frequencies of the order of or larger than H*, and 
low frequency modes with eigen frequencies of the order of H ~ Qp. For the high frequency 
modes, the effects of rotation on the structure of eigen functions and eigen frequencies can 

^ Let us recall that and M* are the radius and the mass of the planet, respectively, and M is the mass of the object exerting 
tides on the planet. 
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be taken into account as a perturbation (see eg. IP) and references therein). Calculation of 
the eigen spectrum of the low frequency modes is much more difficult (Papaloizou & Pringle 
1981, hereafter PP). However, a convenient approximation can also be made in this case 
which allows one to reduce the general equation m to a much simpler equation for the low 
frequency modes only (PP, PI). 


3.1 Eigen problem for low frequency modes 


In order to reduce equation 03 to an expression appropriate for low frequency modes it is 
convenient to use the potential (see equation (ED for definition) as a dynamical variable 
instead of the displacement vector We express the components of in terms of Wm 
using the Fourier transformed variant of equation having set /j, = 0 there. For a normal 
mode labelled by k, the components of the displacement vector are related to the potential 
through: 


^■uu,k 

^z,k — 


(^kdk \ 

^kdk \ 
1 dWk 


al dz 


dWk , 2mH ^ 

-1-hCfc 

ozu w , 

akm 


VJ 


Wk + 2Vt 


dWk 

dw 


(38) 

(39) 

(40) 


where dk = — o’l- To simplify notation we shall not use tildes and subscripts denoting 

m for quantities associated with eigen modes. 

Now we substitute equations (EHD, (Ei and iOD into right hand side of the continuity 
equation Q. The density perturbation on the right hand side is expressed in terms of Wk 
and perturbed gravitational potential with help of equation where we set 4/®** = 0. 
We then obtain 


alAWk - akBWk - CWu = (jldk^i^f^ - Wk), 


and 

A = 

B = 


1 d f d \ 

w dw \^dw J 

2mf2 dp ^ 
w dw'' 


d ( cl \ w?p 
dz \ dz) w'^ 




(41) 


(42) 

(43) 


Note that the operator A can be represented in an invariant form 
A = -V“(pV„). 


(44) 
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As follows from equation (jH)), A is just a modified Laplace operator. 

The operators A, B and C are analogous to the operators B and C introduced above. 
They are self-adjoint with respect to the inner product ©• Also when W is not a constant 
A and B are positive definite and C is non negative dehnite. 

Note that an analogous equation used by PI has a factor — a^) instead of aldk 

on the right hand side. Accordingly, the potential used by PI is different from the 
potential Wk used in this Paper: 

Equation for the perturbed gravitational potential is obtained from equations 0 
and 0 


ATf* + 


AnGp 


^T^GP.jAnt _ 

- 9 —^ k ~ 

CZ C 


Wt 


(46) 


Equations (EH) and (03) form a complete set. They are equivalent to equation m and 
can be used instead of it to find the solutions for the eigen modes. 

Assuming that ~ hi it is easy to see that the terms on the right hand side of equation 
m are much smaller than the terms on the left hand side. In the leading approximation in 
the small parameter they can be neglected and we obtain 


trlAWk - aBWi - CWk = 0. 


(46) 


Equation EH can be used for calculation of the low frequency eigen functions and eigen 
frequencies. Then a formalism exactly equivalent to that given by PI is obtained. 


3.1.1 Equation (El) as a self adjoint problem 

Equation EH can be brought into the standard self adjoint form analogous to equation m 
(PI). For that let us consider two dimensional vector with components 

= aA^/^Wk, = C^/^Wk. (47) 

Now it follows from (EH that it can be written in the form 

akYk = HYk, (48) 


where 




(49) 


Since the off diagonal terms in the matrix (EH are adjoint of each other and the diagonal 
terms are self adjoint it is clear that the operator H has the required self adjoint form. 
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The inner prodnct indnced by the operator H is 
« Yk\Yi »= akai{Wk\AWi) + {Wk\CWi). (50) 


It is obvions that different solutions of equation are orthogonal with respect to this 
product. 

From equations (EH), Ei-cinD and dsni) it follows that the norm corresponding to 
the displacement vector associated with Wk is related to the norm Uk = al{Wk\AWk) + 
{Wk\CWk) induced by the inner product as 

Nk = nk/{aldk). (51) 


The density perturbation determined by a solution of equation (Hi follows from equation 

a 

P = -,(W, - ¥p), (62) 

where the gravitational potential is found from equation (EH with Wk as a source. 

The coefficients Sk determining the energy and angular momentum transfer (see equa¬ 
tions EHD and EED) follow from equations (I2S1) and (EH- We have 

Sk = j wdwdz^iWk - ^kW^\ (53) 

where we take into account that we can use real expressions for the potentials Wk and T®**. 

The expression dH can be brought in a simpler form. For that, let us formally solve 
equation (USD 


Tr' = ^T^GAl%Wk, 
ci 


(54) 


where ^ is the inverse of the operator A* = A -|- It is easy to see that A^ ^ is self 
adjoint. We may, therefore, write the last term on the right hand side of (1^ as 

/ vud'cudz^'^k'^^^^ = 4:7 tG [ vud'cudz^{A~^^Wk)'^'^^^ = — f vudvudz^Wk'^^, (55) 

J cj J cj cj J cj 

where the potential 4/^ is determined from equation analogous to equation (HH but with 
negative of the potential 4/®** as a source 

(56) 


A^i + = 


^^Gp.T.ext 


^2 ^2 


The ’polarisation’ potential 4/^ has a simple physical meaning. If the external potential 
4/®’”* is assumed to be static, the potential 4/^ is induced in the star by the corresponding 
perturbation of density. 

Finally we obtain 
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Sk= j + ^')- (57) 

4 EXPLICIT EXPRESSIONS FOR THE ENERGY AND ANGULAR 
MOMENTUM TRANSFER 

In this Section we obtain explicit expressions for the transfer of energy and angular mo¬ 
mentum in terms of quantities characterising the perturbing tidal field and quantities char¬ 
acterising the planet. These will be applied directly to the problem of evaluating the tidal 
response during a fly by or as a consequence of multiple encounters during circularisation 
from a highly eccentric orbit. Since the expressions for the energy transfer for the funda¬ 
mental and p modes are well known (see PT, PI, IP and references therein) we consider only 
the low frequency modes in this section. 


4.1 Fourier transform of the tidal potential 


As we have mentioned before we consider highly eccentric orbits in this Paper, with eccen¬ 
tricity e ~ 1. In this case the energy and angular momentum transfer takes place mainly 
near the periastron and we make the usual approximation of the orbit as a parabolic one 
with the same value of Vp for calculation of the perturbing potential 

The Fourier’s transform of the perturbing potential for the parabolic orbit in the 
quadrupole approximation has been calculated by PT. We have 


^ m 


5 (2 — \m\)\ B\ 


\m\ 


h-m (cos 6 ), 


\ 27r3 (2 -F \m\)\ (1 -h q) 


(58) 


where B 2 = Bq = q is the ratio of M* to the mass of the star exerting the 

tides, M: q = M*/M. In the application of our formalism to the problem of circularisation 
of the extra solar planets q ~ 10“^ is very small and will be neglected later. 


GM{l + q) 

P 

is a typical frequency of periastron passage. The functions 

POO 

h-m{y)= / da;(l-1-a;^)“^ cos[\/2j/(a;-1-x^/3)-I-2marctan(a;)] (60) 

Jo 

specify the dependence of on a, they are described by PT. y=r2(ct -|- m), where we use 
dimensionless frequencies 




a = cr/fl, n = VL/Vtp. 
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P[^{x) is the associated Legendre function, and (r, 0) are the spherical coordinates. 

When 2/^1 the quantities l 2 -m{y) decay exponentially. This is valid for the case of the 
high frequency fundamental and p modes where 3> 1 for the eigen frequencies. On the 
other hand, the low frequency modes considered above have ctfc ~ 1 and do not experience 
this suppression. 

4.2 First passage problem 

The expressions for the energy and angular momentum transfer ()d4|l and ()d6|l are strictly 
valid only under assumption that the planet is unperturbed before the fly by, and we call 
the problem formulated in such a way as a ’first passage’ problem. This assumption may not 
be correct when the orbit of the planet is elliptical and the viscous time ty = 1/is larger 
than the orbital period Porb- In latter case the energy and angular momentum exchange 
between the orbit and the modes depends on the state of the pulsational modes before the 
periastron passage. However, the expressions for the hrst passage problem provides a basis 
for a treatment of that more complicated case, and hrstly we would like to discuss this 
problem. The generalisation to the case of a multi-passage problem is discussed in the next 
Section. 

In our discussion below we assume that effects determined by perturbations of the struc¬ 
ture of the planet due to rotation are not significant for our purposes. We use, accordingly, 
spherical models of the planet in order to hnd the form of the density and the sound speed 
entering in equations dsoD, (jsni). In the same approximation, taking into account that the 
perturbing potential is proportional to the product of a function of the radius r = 
and a spherical harmonic, we see that the variables in equation (jini) can be separated in 
spherical polar coordinates. Therefore we can express the sum entering in equation 

flhTjl as 

^,ext + (62) 

where the function F[r) can be obtained from equation (jSSI)- 

In order to present the expressions dsoD, dsni) and dSZl) in a simpler form we express the 
quantities Ej and Lc in terms of natural units of energy F* = GM^/R^ and angular 
momentum {GM^R^). As above we also adopt dimensionless units for the spatial 

coordinates, density and the sound speed expressing them in units of i?*, the averaged 
density p = and respectively. 
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Instead of the separation distance Vp it is convenient to introduce the parameter (PT) 


T] = = 


M* r3 


\{l + q)M R^- 

We then substitute equations and into (On|l and (pn|l . We obtain 


F 

™ (l + g)2 


k 


Lo — 


2C. 


2 o3 


- 4)Qlkl-Jyt)y: 


(63) 


(64) 


(65) 


(1 + g)" .. ?7^ 

where we separate contributions to Ec and Lc corresponding to different values of m such 
as Ec = J2m=o,2 Em and take into account that Lq = 0 and therefore Lc = L 2 . All quantities 
in equations dnn), (ESI) are assumed to be in the dimensionless form, (72 = ^, Co = 3/4 and 

= fl(d-(fc) + 2). (66) 

The overlap integrals Qk are analogous to the same quantities introduced by PT in 
their analysis of dynamic tides in a non rotating star. They describe spatial coupling of the 
perturbing potential and of the eigen modes, and have the form: 

Qk = {^E{ryP^\Wk)/V^k, (67) 

where fbk = nk/Ll^. 

The expression (El and (El have already been obtained by PI using a different method. 
They can also be presented in another equivalent form as 


Cm(f^) 


(1 + g) V 


77* 


Lo — 


A(0) 


(l + g)V"*’ (68) 

where the dependence of the functions and A on hi follows from El and (El- Accordingly, 
the energy transfer in the inertial frame follows from El, El and the use of 


Cm = - ^lk))Qlyl-m{yk), 

(k) 

where all quantities in the brackets belong to a particular value of m = 0 , 2 , 
A = 2C2ri^ ^ ct^;i.)(4 — d'‘^k))Qlk)^2,-2iyk)^ 

(k) 

where all quantities in the brackets are evaluated for m = 2 and 
6/ = 62 + hlA, 


as 


Er = 


(eo + 62 + A) 

(1 + g)2?7® 


77* = 


(ep + ej) 

(1 + g)2?7® 


77* 


(69) 


(70) 


(71) 


(72) 
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Usually, the moment of inertia of the planet is much smaller than the moment of inertia 
of the orbit. In this case the planet tends to rotate with ’equilibrium’ value of the angu¬ 
lar velocity, Qeq determined by condition that there is no exchange of angular momentum 
between the planet and the orbit at this rotation rate, 

A(fleg) = 0. (73) 

This state of rotation is often referred to as a state of ’pseudo-synchronisation’. For a given 
model of planet, equation dZSD can be solved numerically to hnd fleg- The expression for the 
energy transfer at the state of ’pseudo-synchronisation’, Epg, follows from equation (IZ2D 

Eps = ^E,, (74) 

where e* = depends only on the structure of a planet. For simplicity we set 

1 -|- g = 1 hereafter. The expressions for a comparable mass ratio q can be easily obtained 
from our equations by an appropriate re-scaling of variables we use. 

An estimate of e* can be obtained from the following semi-analytical arguments. As we 
will see below only two ’global’ modes with m = 2 mainly determine the energy transfer 
when ~ fleg and we can take into account only these two modes in the summation cza) 
and (IZID- One of these modes has a positive eigen frequency with ui ~ 0.6 and the overlap 
integral Qi ~ 0.16, and another mode has a negative eigen frequency with ct 2 ~ —1.1. 
The equilibrium value of the angular velocity is always close to the frequency of periastron 
passage, with Oeg ~ 1.55. This information is sufficient to estimate e*. Indeed, we can express 
the second overlap integral from equation (na in terms of Qi, the eigen frequencies and the 
quantities /2,-2(2/1,2), where 1/1^2 = ^eq{‘^ + (^i, 2 )- We get Q 2 ~ 0.033. Taking into account that 
^2,-2(2/1) ~ 0.32 and h,-2(1/2) ~ 0.74, we get e* ~ 4 ■ 10“^. This value is in a good agreement 
with numerical results discussed below, see discussion and equation (IHHl) in Section 6. Note 
that for n = 1.5 polytrope e* ~ 6.5 ■ 10“^ (PI)- 

4.3 Multi-passage problem 

In the circularisation problem expressions (Ea, dsa, dsa, dni, (ina, dna, dza and jza 
are strictly valid only when the energy stored in the modes can be dissipated between two 
successive periastron passages. In the opposite limit, when the dissipation time tu 3> Porb, 
where Porb is the orbital period, the planet approaches periastron in a perturbed state. 
The additional linear perturbation induced during a particular periastron passage must be 
added to the perturbation already present in the star. Since the mode energy is quadratic 
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in the perturbed quantities, it can either decrease or increase after the periastron passage, 
depending on phase of the perturbations (Kochanek, 1992, Mardling 1995a, 1995b). This 
issue has been studied analytically by IP. It was found that when a change of the orbital 
period between two successive periastron passages is sufficiently large, the energy transfer 
proceeds in a stochastic way. In this case the phase of perturbation in the planet approach¬ 
ing periastron is not correlated with the phase of the perturbation induced by tides near 
periastron, and the energy transfer during the periastron passage is essentially a random 
quantity, with dispersion of the order of the expressions for the energy transfer calculated 
above. In this situation, a build-up of mode energy due to stochastic instability is possible. 
A semi-analytical condition for this instability to happen has been obtained in IP. Here we 
generalise their formalism to take account of the specihc features associated with inertial 
waves. 

Let us assume that dissipation is negligible and treat the tidal influence as a sequence 
of impulses separated by periods of time equal to the orbital period. In this case, any time 
between two impulses the planet undergoes free oscillation, and the Lagrangian displacement 
vector can be represented by an expression analogous to (EZD 


^ = 27r2 ^ 


CTkSk 






(75) 


m,k 

where the dimensionless real coefficients ^ 1 characterise the amplitudes of perturbations 
corresponding to different modes and ijjk are the corresponding phases. Note that the origin 
of time is chosen to coincide with the moment of periastron passage. The expressions for the 
canonical energy and angular momentum follow from (EHl), (|S3) and m 


E, = 8Tr^J2^ 


2^lSkSl 


Nk 


Lc = 167r^ 


SkSl 

Nk ■ 


(76) 


m,k m=2,k k 

When Afc = 1 these expressions obviously coincide with the expressions (P|) and for 
the single passage problem. For inertial waves equations (I7U|) can be brought to the form 
analogous to equations (IHHl-irnil 






Ej, 


(1 -|- ’ 


-^2 — ( X! ^k^l 


(1 -I- qY'rj^'’ 


(77) 


where we now show separately the contributions of terms with different m to the canonical 
energy. em,k and Xk are the terms in the summation series for Cm and A, corresponding to 
a particular value of k, see equations dHED, (innD and (HOI). As we discus below in the tidal 
problem we can neglect the contribution of m = 0 terms. In this case, the mode energy in 
the inertial frame can be represented as 
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E, = I ^ 




ps? 


(78) 


\.m=2,k 

where Epg is given by (I7il) 

After the periastron passage, the expression (EZD must be added to m, and the resulting 
expression can brought again to the form dZSD but with new values of the amplitudes and 
phases, and -0^. Thus, the action of tidal forcing dehnes an iterative map —»• 

(A^,'0fc)- It is convenient to introduce new variables cos cos V’fc) 

and {yk~^,yk) = (A^“^ sin-^^”^, A^sin-^^). In terms of these variables this map has espe¬ 
cially simple form (IP) 

/ 


Xi 


= R{-(t)k) 


+ 1 


Vk ^ 


(79) 


V Vk 

where (f)k = (^kPorbi and -R(0) is a two dimensional rotational matrix. Note that in equation 
CSl) we take into account the change of origin of time after the periastron passage: t —>■ 
t “1“ Porb' 

When Porb is a constant, it is easy to see from equation (HSl) that the amplitudes A^ do 
not grow with the number of iterations. However, the orbital period depends on the orbital 
energy, Eorb, as 


Porb 


(80) 


\J‘^\Eorb\ \Eorb\{^ + 

and the orbital energy is proportional to the mode energy Ej due to the law of energy 
conservation. Since Ei = ^ps is proportional to the squares of the 

modes amplitudes, A^“^, the map is non-linear. 

In order to simplify the problem, let us consider a situation when the mode energy is 
much smaller than the orbital energy such that the orbital energy and angular momentum 
remain in the neighbourhood of some particular hxed values. In this case we can set Porb = 
-Porb(O) -l- APorb, where Porbi^) is a constant, and 

a5 El 


AD ^ P 

/XiQrb —--r , ^ — —OTT 

2 Eorb M, 


\ {GMfMY 


(81) 


where Eorb and a are the orbital energy and semi-major axis corresponding to the period 
Porb{0). Now the angles <pk can be expressed in terms of Ej as 


4>k = 4>k{0) — akCx-Ej, 

where Ei = Ej/Eps, 0^(0) = akPorbiSY and 


(82) 


© 2002 RAS, MNRAS OOP.ITICT 











Dynamic tides in rotating objects 21 



Figure 1. The dependence of acrit on 0i(O). The solid and dashed lines correspond to the maximal number of iterations 
Nmax = 10^ and Nmax = 10^, respectively. 


u] ICD 



lOXW 


Figure 2. Ej as a function of the number of iterations N. The case (/>i(0) = 2.32 corresponding to the maximal value of 
acrit = 1-04 is shown. The solid curve corresponds to q = 0.9 and the dashed curve corresponds to o; = 1.1. 


a = 67re*0 


M* 1 / 


(83) 


KtTrJ' 

where we use equations (ED), (El, dZl and m- Equation (IHl is analogous to equation 
(84) of IP. 

Equations dlHl), dH and (El fully determine the map for given values of 0^(0) mod 2 tt 
and a. When a are sufficiently large, stochastic instability sets in. Since the angles 0fc(O) 
are determined by initial conditions, the results should be averaged over these angles. 

This map has been investigated numerically by IP for the case of only one high frequency 
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mode. It has been found that the instability sets in when the average value of a is larger than 
ttcrit ~ 0.67. In our case the situation is slightly more complicated since there are at least two 
modes which play a role in the energy transfer. Also, we should specify the evolution of the 
angular velocity. In order to investigate this question we consider numerically the simplest 
possible case of constant equilibrium rotation of a planet, when (2 = For dehniteness, we 
use parameters described in the end of the previous section, namely fleq = 1.55 and consider 
only m = 2 modes with cti = 0.6 and d 2 = —1.1. In this case we have 62 ~ 0.35e* in the hrst 
case and 62 ~ 0.65e* in the second. We iterate numerically the map for a large number N of 
iterations for a given value of 0i(O) with 02 ( 0 ) = |^0i(O), and for a given value of a. The 
stochastic instability is assumed to set in when the energy Ej exceeds the maximal number 
of iterations, N^ax during the iterations. This happens when a exceeds some critical value 
'^crit(0i(O)). The results of our computations are shown in Figure 1. One can see from this 
Figure that there is a strong dependence of acrit on 0i(O). Nevertheless we use the averaged 
value of acrit: ^crit = ^ So"" otcrit ~ 0.43 since our results are rather insensitive to the exact 
values of acrit- It is interesting to note that the transition from a quasi-periodic behaviour 
of E] to stochastic behaviour is very sharp. In Figure 2 we show the evolution of Ej for 
01 (0) = 2.32 and two values of a: a = 0.9 < acrit = 1-04 and a = 1.1 > acrit = 1.04. One 
can see from this Figure that a small change of a is sufficient for stochastic instability to set 
in. 


From the condition a > acrit we hnd that the stochastic instability sets in only when 
the semi-major axis is sufficiently large. 


a ^ Ost ( 




67rf2p^6^ 


\Mj 

Using the values of all quantities in 
solar planets, we have 


(84) 

typical for the problem of circularisation of extra 


^st 


30 hrr 

VM* 


R 


Mr, 




14/5 


au. 


(85) 


where rjiQ = rj/10, Mj ^ 2 ■ 10^°and Mq ^ 2 ■ 10^^ are the masses of Jupiter and of 
the Sun, respectively, and Rj ^ 7 ■ lO^cm is the radius of Jupiter. Remarkably, the typical 
value of Ost ~ 30aM is close to what has been obtained in the analysis of the stochastic 
instability associated with the fundamental modes (IP). This value is of the order of a 
typical ’initial’ semi-major axis from which we would expect the process of circularisation 
to begin for rj = 10. This corresponds to a hnal circularised period of 3.3 d for a Jupiter 
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mass planet around a solar type star. However if the final period is reduced to ~ 1,2 d, 
corresponding to the shortest period exoplanets, the initial semi-major axis is reduced to 
between 1 and 2au. If the initial semi-major axis is smaller, the stochastic growth of the 
mode energy is not possible. However, in this case, we would expect a stronger dissipation 
of the mode energy due non-linear mode-mode interactions (e.g. Kumar & Goodman 1996). 
This possibility needs a further investigation. 

5 PLANET MODELS AND THEIR NORMAL MODE SPECTRA 

In this Section we construct realistic models of planets at different evolutionary stages for 
which we obtain the oscillation spectra numerically. We discuss our numerical method, use 
it to solve equation ii, and discuss the results which are used in subsequent sections to 
evaluate the tidal responses of different planet models undergoing a fly by or circularisation 
from a highly eccentric orbit. 

We hnd it convenient to adopt a system of dimensionless units for the spatial coordinates, 
density and the sound speed expressing them in units of i?*, the averaged density p = 
and respectively. We also use dimensionless frequencies 

ct = (t/H, H = H/Hp. (86) 

5.1 The planet models 

We assume that the planets are spherically symmetric and adopt the equation of state of 
Saumon, Chabrier & Van Horn (1995). Although it is likely that giant planets form through a 
nucleated core instability (eg. Bodenheimer & Pollack 1986), for this hrst set of calculations, 
for reasons of simplicity we shall assume that solid cores are absent. In this context we note 
that the mass of any solid core in Jupiter is highly uncertain, being sensitive among other 
things to the equation of state adopted, with recent estimates showing it to be consistent 
with zero (see Saumon & Guillot 2004). To make this consistent with a nucleated core 
instability scenario for formation, these authors suggest that a core could be eroded and 
remixed during subsequent gas accretion of light elements. But it should be noted that the 
presence of a large solid core could signihcantly affect the mode spectrum of the planet 
through the existence of wave attractors (eg. Ogilvie & Lin 2004) that could produce effects 
on a long time scale. However, the presence of a small core is unlikely to affect phenomena 
of short duration such as the energy transfer that occurs during a close passage or fly by. 
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Figure 3. The density p scaled in terms of the mean density as a function of the dimensionless radius r. The three solid curves 
are for models with M* = IMj and radii R* = IRj, lARj and 2Rj respectively. The three dashed curved are for models with 
M* = 5Mj and radii R^ = lARj and 2Rj respectively. The dotted curve corresponds to the model with M* = IMj 

and radius il* = IRj. and ’smoothed’ density profile. Curves of the same type starting from larger values at r = 0 correspond 
to planet models with larger radii. Note that because of the use of different scalings to dimensionless variables the numerical 
values differ from those of IP by a factor ^. 



Figure 4. As in Fig. |3 but the absolute value of the density gradient |^| is shown. Curves of the same type starting from 
larger values at r = 0 correspond to planet models with larger radii. 

We calculate numerically the form of p and Cs as functions of the radius r. With these 
functions known, equation (I48j) can be solved and the overlap integrals (EZD calculated. As 
in IP we consider only two planet masses M* = IMj and 5Mj. 

For a given mass we calculate planet models for six values of the radius. These are 

= IRj (i?* = 1.03i?j) when M* = Mj (M* = 5Mj), and i?* = 1.2Rj, lARj, 1.6Rj, 
1.8Rj, 2Rj for both masses. The different values of radius for a given mass correspond to 
different points on an evolutionary track. On such a track, the radius decreases with time 
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Figure 5. As for Fig. and Fig. but the square of sound speed is shown. Curves of the same type starting from larger 
values at r = 0 correspond to planet models with larger radii. 

due to gravitational contraction. To determine the age corresponding to a given radius, 
we use the relationship between the radius and luminosity, and hence the age of the planet 
calculated by Burrows et al (1997) (see Fig. 8 of IP). Since a young hot planet of age ~ Itfyr 
has a radius of around ~ 2i?j, the range of radii we consider corresponds to planet ages 
between ~ Itf'yr and ~ 5 ■ 10® — lt)^^yr. 

We illustrate the dependences of the density, the magnitude of the density gradient | ^ | 
and square of the sound speed for several planet models in Figures OE As one can see from 
Figure El for a particular planet, the ratio of central density to mean density increases with 
radius. This effect is less prominent for the models with M* = 5Mj. 

Note a sharp drop of density ai r k, 0.82 when M* = IMj, and = IRj. This is due to the 
phase transition of hydrogen from the metallic to molecular phase which is assumed to be 
of the hrst kind. This phase transition also takes place when M* = 5Mj and i?* = 1.03i?j, 
but it is shifted to a larger r 0.96 and is not well resolved in this Figure. In general, 
according to Saumon, Chabrier & van Horn (1995), the phase transition is of the hrst kind 
when i?* < 1.27i?j (i?* < 1.255i?j) for M* = IMj (M* = 5Mj). For larger radii the adiabats 
describing the interiors of the planets go above the critical point of Saumon, Chabrier Sz 
van Horn (1995) equation of state on {P,T) phase plane. As we see below the effect of the 
phase transition on the structure of the eigen spectrum is more prominent for M* = IMj. 
To illustrate this it is of interest to compare this model with a similar model without the 
phase transition and we make such a model with ’smoothed out’ transition between the 
metallic and molecular phases (dotted curves in Figures 3-5), for M* = IMj and R^ = IRj. 
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This model has been calculated using a prescription provided by Saumon, Chabrier & van 
Horn 1995. 

The structure of the density profiles is better seen in Figure E] where the magnitudes of 
the density gradients are plotted. For the models with the phase transition there is a dis¬ 
continuity at the radius of the phase transition. For the models without the phase transition 
the modulus of the gradients increase with radius at small r, and they sharply decrease at 
large r. The transition point between these two limiting cases corresponds to a ’melting’ 
point where the number density of the hydrogen molecules sharply decreases toward smaller 
values of r. This transition point shifts toward smaller r with increase of the planet radius. 
As seen from Figure El the sound speed is changes non monotonically near the transition 
point. 

In general, the models corresponding to the larger mass are more similar to each other. In 
this case the radius of transition is shifted toward larger r where the planets have rather low 
density. The low density outer regions of the planet should not influence eigenfrequencies and 
eigenfunctions corresponding to global modes should they exist. Therefore, it is reasonable 
to suppose the presence of the phase transitions and the ’melting’ zones in the planet models 
with = 5Mj is less significant for our problem than for the case of M* = IMj. As we 
see below this assumption is supported by our numerical results. 

5.2 Numerical method 

For our numerical work we use a method similar to that proposed by PP. This involves 
reducing the operator equations to discrete form by use of a set of basis functions. 

5.2.1 Matrix representation of equations and 

For our numerical work it is convenient to work with a matrix representation of equation 
(jini). This is accomplished by use of a basis of ’trial’ functions v\ j = 1, ..cxo which is assumed 
to be complete. We introduce a Fourier decomposition of an eigenfunction W(^k) in terms of 
the trial functions as 

= (87) 

3 

where the indices enumerating the eigen functions are enclosed in brackets from now on. 
Then we substitute equation into equation (HH|l and taking the scalar product Q of 
the resulting expression with n* we get 
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( 88 ) 


- C|)af' = 

where 

A{ = {D\Av^), Bi = {D\Bv^), Cl = {v^\CC), (89) 


and the standard summation rule over the repeating indices is implied We can introduce 
square roots and inverses of the above matrices in a similar manner to that adopted for 
the corresponding operators. Let us consider inhnite dimensional vectors analogous to the 
vectors F(fc) (see (jUj)), with components where the index M stands for the pair ((l,i) 
and (2,i)). We have 


y,'‘> = y:W = (90) 

The eigen value problem ()88j) can be formulated for the vectors o in a way equivalent to 

(gSt, 

a,t,yy> = (91) 


where the matrix has a structure equivalent to the operator H (equation ()49|l l with 
the operators being substituted by the corresponding matrices. 

It is natural to use the solutions of the eigenvalue problem 

Av^ = Xpv^ (92) 


as the functions vA These functions are orthogonal with respect to the inner product 0. 
As follows from equation (I44|l . for the spherical planet models the variables in equation 
are separable in spherical polar coordinates {r,6), = Vn,i{r)Pl^{cos 6), where n = 1,2,.... 

n — 1 is the order of the radial eigenfunction and the index j stands for the pair of positive 
integers (n, /). For the radial part we have 


WW + ^P{r)Vn,l = Xn,ip{r)vn,l, 




dr 


(93) 


where L = ^Jl{l + 1). It follows from that the eigen values of A are bound from below: 
Xn,i ^ when I ^ 0. 


5.2.2 Numerical approach 

In our numerical work we truncate the series (IHID at some sufficiently large value j = jmax- 
The trial functions must be arranged in such a way that when j is sufficiently large the 

^ We stress that we do not sum over the indices enclosed in the brackets that label an eigenfunction. 
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corresponding trial functions oscillate with periods of spatial oscillation that decrease with 
increasing j. In this case, one can argue that errors will be small because of a negligible 
contribution of the oscillatory terms with high j to the various integrals occurring in the 
formulation of our problem as discussed above. Of course, this idea must be tested, and we 
discuss such tests later. 

Since for our problem m = 0, 2 is even and the angular parts of the trial functions are 
such that Pi^{x) = {—iyPl^{—x), the modes odd with respect to reflection in x —>• —x do 
not couple with the tidal held and we can choose Z ^ m to be even. Thus I = 2li with Zi 
being a positive integer such that li ^ m/2. 

It is convenient to consider an equal number of the radial functions n{max) and angular 
functions li{max). Thus we write n{max) = li{max) = Nmax- In our numerical work, we use 
different values of N^ax in order to study the dependence on the number of trial functions, 
with the largest N^ax = 15. Accordingly, we have jmax = N/^ax ^ 225. 

After truncation of the basis set used in (IHZD we have ordinary hnite dimensional matrices 
Al, Bl and C/. We calculate the eigenvalues eigenvectors of these matrices. The eigenvalues 
are used to calculate the inverses and square roots of the matrices entering in the matrix . 
We hnd the components of and solve the eigenvalue problem (jH^. The eigenvectors 
are used to hnd the coefficients af'\ and we use these coefficients to hnd the eigenfunctions 
W(k) with help of equation (jHZl). 

One can view the eigenvalues and eigenfunctions as belonging to the oscillation spectrum 
of the planet. For the case of a model without a solid core this is ultimately likely to 
be discrete but everywhere dense (see eg. PP). One should also consider the possibility 
of the existence of a continuous spectrum in some cases. Then the numerically obtained 
eigenfunctions may not show convergence to regular functions. But note that, independently 
of this issue, we expect our numerical method to be convergent when the eigenfunctions we 
obtain are used to provide a basis for representing smooth functions varying on a global scale. 
The issue of the possibility of ultimately irregular eigenfunctions belonging to a continuous 
spectrum does not ahect this or preclude the convergence of the process of calculating the 
angular momentum and energy exchanges associated with a fly by. On account of the forcing 
by a global potential and short time duration of a fly by, a basis representation capable of 
representing global scales would be expected to be adequate on physical grounds. Indeed for 
our models it turns out that only a few global eigenfunctions are important for determining 
the energy exchanged and these can be represented with a basis set that needs only to 
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be capable of representing global fnnctions. This makes the issne qnalitatively different to 
the one that occnrs in the context of tidal forcing occnrring in a circular orbit with hxed 
frequency for inhnite time periods (eg.Ogilvie & Lin 2004). 

It is obvious that the dimension of the matrix is 2jmax- Thus, we have 2jmax eigen 
values (T(fc) and eigen vectors The number of terms in the series occurring in (IMI) . (10511 and 
dZZI) determining the tidal response, is proportional to jmax- If is important to check whether 
the process adopted and the corresponding series converges and we discuss this issue below. 
As we have stressed above, the convergence of the series means that the tidal response is 
mainly determined by a few ’global’ modes with a large scale distribution of W{k) over the 
planet. 

5.2.3 Degenerate modes with a(^k) = 0 and the case with m = 0 

The results of our numerical calculations reported below allow us to suggest that the nu¬ 
merical method described above is stable and robust for the modes with sufficiently large 
values of |o'(fc)|. However, there is a problem with this approach related to the case of small 
values of |cr(fc)| 0. This problem is associated with the existence of a non-trivial null space 
of the operator C. It is easy to see that all modes satisfying 

CW^k) = 0 (94) 

are degenerate solutions of (P6|l (and with cr(fc) = 0. In fact any W(k) that is independent 
of ^ will do. The corresponding norm n(^k) = 0. Physically, for m = 0, these modes correspond 
to introducing an inhnitesimal amount of differential rotation and thus slightly shifting the 
static equilibrium. As follows from equation (EZD the overlap integrals are formally inhnite 
for these modes. However, in the exact theory this does not lead to any difficulty since the 
contribution of a particular mode to the energy exchange is proportional to (see e.g. 
equation This means that the degenerate set of modes is in fact not excited. 

However, due to numerical inaccuracy of our method the modes formally belonging to 
the null space of C acquire small non zero a(^k) ■ The corresponding overlap integrals are very 
large and they can lead, in principal, to a very large unphysical energy transfer. This problem 
is not signihcant for the case of m = 2. In this case, the quantities l 2 ,- 2 {y) entering in (15411 
are very small for an interesting range of y. They suppress effectively the contribution of the 
spurious modes to the energy and angular momentum exchange. On the other hand, in the 

^ Note that some of <y[k) degenerate. 
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Figure 6. The overlap integrals Q(fc) as functions of the eigen frequency The dotted and dashed curves correspond to 

jjnax = 25 and jmax = 100 respectively. Circles and squares show the positions of particular Q(^k)- The solid curve corresponds 
to the case with jmax = 225. 

case of m = 0 the corresponding qnantities hpiv) do not snppress this contribntion, and the 
spnrions modes do lead to a very large nnphysical energy transfer. Therefore, in principle, 
the nnmerical method must be able to overcome this difficulty, by e.g. incorporating the 
solutions to (ED in the set of the trial functions. 

This can be done by considering polynomials consisting of products of different power 
of zu and 2 ; as a new set of trial functions (PP) instead of solutions to El . However, this 
set is difficult to use when considering a large values of jmax- Therefore, in our numerical 
work we simply identified the spurious modes corresponding to m = 0 by comparing the 
spectra calculated using these two different sets of trial functions and removed them from 
the summation in the expression for the transfer of energy (there is no angular momentum 
exchange in this case). The resulting contribution to the energy exchange appears to be very 
small compared with what is obtained for the case of m = 2, and we will not discuss further 
the case m = 0 below. 

5.3 The spectrum of a planet with M* = IMj and R = IRj 

In this Section we discuss in detail the case of a planet with M* = IMj and i?* = IRj. 
This case has the most complicated spectrum of all models we consider. In our opinion, 
this is mainly because of the presence of the phase transition at a sufficiently small radius 
r 0.82. The analysis of the eigenfrequencies and the eigenmodes developed for this case 
can be easily applied to planets with larger radii and masses which have a simpler structure. 
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Figure 7. The quantities ej, €2 and A as functions of the dimensionless rotation rate Q. The dotted, dashed and solid curves 
correspond to jmax = 25, 100 and 225, respectively. The curves which always pass through positive values represent e 2 . The 
curve of a given type that takes on the least value corresponds to ej, the remaining curve corresponding to A. 





Figure 8. The binned overlap integrals Qb,(k) plotted for jmax = 25 using the dotted curve with circles; for jmax = 100 
using the dashed curve with squares; and for jmax = 225 using the dot-dashed curve with triangles. The solid curve shows the 
original overlap integrals Qn^(^k) calculated for the case of jmax = 225. Note that the dot-dashed and the solid curves almost 
coincide. 


In Figure IHl we show the overlap integrals Q(fc) (see equation (IF)7|l for the dehnition) 
calculated for this case as a function of the eigenfrequency d'(fc). It can be shown analytically 
(PP) that |ct(fc)| < 2. The curves of different types correspond to different numbers of the 
trial functions, jmax = 25, 100 and 225. We show positions of particular Q(^k) representing 
results of calculations with jmax = 25 and jmax = 100 as circles and squares, respectively. As 
we have mentioned above the sharp rise near cr(fc) ~ 0 is due to numerical inaccuracy. There 
is a similar unphysical rise near a(^k) ~ 2. It is clearly seen from Figure IHl that the width 

© 2002 RAS, MNRAS 000. [TII^ 


























32 P. B. Ivanov and J. C. B. Papaloizou 



Figure 9. Same as Fig. |7|but now the dotted, dashed and dot-dashed curves are calculated using the binned overlap integrals 
for jmax = 25, jmax = 100 and jmax = 225, respectively. The solid curve is calculated using the original overlap integrals for 
the case of jmax = 225. 



Figure 10. The overlap integrals Qn,k ^re shown for two models of the planet: the model without the phase transition - dashed 
curve, and the original model with the phase transition - solid curve. 


of these numerical features depends on the number of trial functions jmax decreasing with 
increase of jmax- Similar to the case of a polytropic star with index 1.5 (PI), there are two 
’global’ modes with large values of the overlap integrals, with ct(i) ~ 0.5 and ct( 2 ) ~ —1.1. 
The frequencies and amplitudes of these modes depends only weakly on the number of trial 
functions, and these modes are especially important for the problem on hand. However, 
contrary to the case of a polytropic star with index 1.5, there are some other markedly large 
values of Q(^k) associated with modes at other eigenfrequencies (eg. near ~ —0.9 and 
(J 4 ~ 1.25). These can also be associated with global modes, which are, however, less stable 
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Figure 11. Same as Figures f71 and IHI but now the dotted and dashed curves are calculated with using the binned overlap 
integrals corresponding to the global modes. We take into account only the two ’main’ global modes to calculate the dotted 
curve and the eight global modes discussed in the text to calculate the dashed curve. The dot-dashed curve is calculated using 
the overlap integrals corresponding to the model without the phase transition. Note the dotted and dot-dashed curves for ej 
and A practically coincide. 



Figure 12. The distribution of over the planet. W(^k) = 0 at the origin. Contour lines of different style correspond to 

values of of opposite signs (which are always found to occur). Ten contour levels of Wn (n=l,10) for each sign are defined 
as follows: Wn = “fF*, where W* is one of the maximum or minimum value of The upper left (right) plot corresponds 

to the main global mode Wi with cri = 0.521 {W 2 with a 2 = —1.11). The lower left and right plot show the same distributions 
but calculated for the planet model without the phase transition. 
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Figure 13. Same as Fig. II 2l hnt for Ws with 0-3 = 0.873 ( upper left plot); W 4 with 0-4 = 1.257 ( upper right plot); W 5 with 
<75 = 1.273 ( lower left plot); Wq with ae = —1.525 ( lower right plot). 


with respect to changes of jmax- They cannot be reliably identihed with a small nnmber of 
the trial fnnctions {jmax = 25 in onr case). We discnss properties of these modes later in 
this section. 

In Fignre [7| we show the qnantities 62 , A and ej dehned by eqnations (159117111 which 
determine the energy and angnlar momentnm transfer throngh eqnations ()58jl and (1Z2D. 
As for fignre IHl the cnrves of different types correspond to different valnes of jmax- It is 
seen from this Fignre that when O is snfficiently large the cnrves representing large valnes of 
jmax = 100 and 225 are in a very good agreement with each other. The cnrves representing ej 
and A for jmax = 25 are also in qnalitative agreement with those corresponding to the larger 
valne of jmax- However, the cnrve representing 62 deviates signihcantly when Q > 7. When 
H is small, the cnrves representing jmax = 25 and 100 have larger valnes compared to the 
case with jmax = 225. This is dne the contribntion of spnrions modes with eigenfreqnencies 
very close to a = 2. In principle, we can improve the agreement between different cnrves 
signihcantly simply by discarding these modes. However, as we will see below it is possible to 
improve the stability of onr nnmerical scheme with respect to the choice of different nnmbers 
of trial fnnctions and snppress the inhnence of these spnrions modes antomatically. 
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Figure 14. As in Fig. 1121 hut here we show the distributions of W 7 with 0-7 = —1.5 in the upper left plot and Ws with 
(T 8 = —1.353 in the upper right plot. Two non-global modes showing small scale structure with ag = 0.49 (lower left plot) and 
(Jio = —0.5 (lower right plot), are also shown. 


5 . 3.1 Binned eigenfrequencies and overlap integrals 

In order to improve the stability of our numerical scheme with respect to numerical noise we 
would like to separate the eigenspectrum into bins with a hnite frequency width and adopt 
the sum of the overlap integrals within a particular bin instead of Q(k)- To begin, we redehne 
the overlap integral according to the rule 

Qn,(k) = ^(fc) \/(^ ~ ^(k) )Q{k)- (95) 

It follows from equation (IHnj) . that the quantity 62 has especially simple form when expressed 
in terms Qn,{k}, namely 

(2 = C,n*J2Qlf,/i_Jy,). (96) 

( 1 ) 

To dehne the frequency bins let us separate the eigenspectrum into regions between succes¬ 
sive minima of treat all Qn,{k) and a{k) within such a region as belonging to the 

same bin. Note that modes dehning minima are not included in the bins. Accordingly, we 
have 
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- VZ^Wn,(fe)’ ^Ml) - -7^2-’ 

^ Vfe,(Z) 


( 97 ) 


as new ’binned’ overlap integrals and ’binned’ eigenfrequencies. Here, the summation is 
performed over all eigenfrequencies belonging to the same bin, and I is a number of a bin 

The result of comparison between the ’binned’ quantities Qb^k corresponding to different 
jmax and the quantities Qn^k calculated with the largest jmax = 225 is shown in Figure 
IHl One can see from this Figure that the binned quantities calculated for jmax = 100 and 
jmax = 225 are in a good agreement. The results corresponding to jmax = 25 agree with 
what is shown for larger values of jmax for the two ’main’ global modes. They are also in 
qualitative agreement for cT(fc) < 0 , but deviate in a part of the spectrum with positive 
values of a(^k) ■ It turns out that this deviation is not important for our purposes. Indeed, the 
positive part of the spectrum plays a role in the tidal response problem when (2 is small. 
In this case width of the function l 2 - 2 {y) is rather large, and the tidal response is mainly 
determined by the ’main’ global mode with cti ~ 0.5 and the largest overlap integral. When 
0 is large, the tidal response is mainly determined by the modes with negative cT(fc) which 
can then comove with the tidal perturbation near periastron. The width of l 2 - 2 {y) decreases 
with increase of ( 2 , and, when the rotation rate is large, the tidal response is determined by 
several modes with sufficiently large overlap integral and sufficiently large absolute values of 
cr(fc). This effect explains in part why the energy transfer in the rotating frame determined 
by the quantity €2 tends to some non-zero value with increase of the rotation rate (see e.g. 
Figure 7 and Figure 9). 

In Figure El we show the quantities ej, 62 and A calculated using the binned overlap 
integrals for the different values of jmax and the corresponding quantities calculated using 
the non-binned overlap integrals with jmax = 225. Comparing Figure[7|with Figure|Hlwe see 
that the use of the binned quantities does allows us to improve signihcantly the stability of 
our numerical scheme with respect to different numbers of basis functions. Now, all curves 
are in a good agreement with each other for the whole considered range of f 2 . 


5.3.2 Eigenspectrum of the model without the phase transition 

In Figure^! we show the result of calculation of the overlap integrals for the model with 
’smoothed out’ phase transition (dashed curve) together with the overlap integrals for the 

^ Note that this procedure is not unique, and better results may be achieved with a more sophisticated scheme of binning. 
However, our procedure is, perhaps, the simplest one, and it is definitely sufficient for our purposes. 
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model with the phase transition. As it is seen from this Fignre the two ’main’ global modes 
have practically the same overlap integrals for the both models. However, other prominent 
spikes present in the model with a phase transition are absent for the model withont the 
phase transition. Therefore, there only two global modes in that case. One can assnme 
that the presence of the other global ’non standard’ modes in the realistic planet models is 
related to the presence of the phase transition. In fact, as we will see below these modes 
are also present in the planet models with larger radii where the phase transition is absent. 
It is reasonable to snppose that the presence of some of these modes is related to a sharp 
decrease of the density after the radins of the phase transition (or the melting point for the 
models with larger radii, see FignresElandlD) in the realistic planet models. 


5.3.3 Global modes 

As we have pointed ont above there is a set of modes characterised by relatively large 
valnes of the overlap integrals. These modes are referred to as global modes. Physically, 
these modes have eigenfnnctions for which W{^k) varies on a large scale. In this sitnation, the 
overlap integrals are not averaged to small valnes after integration over the volnme of the 
planet. 

As can be seen from Fignres El IHl and uni there are more than ten modes with markedly 
large overlap integrals. However, they are not eqnally important for the problem on hand. 
The two most important modes are ’the main’ global modes with the largest valnes of the 
overlap integrals. These are present in all onr models as well as in polytropic stars. Other 
global modes with negative valnes of (j(^k) can also play an important role in the tidal response. 
Additionally, it is interesting to consider the positive and negative freqnency global modes 
with the largest valnes of Qn,{k) after the main modes. These modes contain information 
abont the difference of the realistic planet models from over-simplified polytropes, and may 
be of interest for another problems as well. There are eight modes satisfying these criteria. 
Their eigen freqnencies and overlap integrals are shown in Table H Wi and W 2 are the two 
’main’ global modes. W 4 and IV 5 as well as Wq and IV 7 have very close eigen freqnencies. 
Most probably, these modes are donblets split as a resnlt of a pertnrbation. 

In Fignre 11 we show the qnantities e/, €2 and A caicniated with only two global modes 
(the dotted cnrves), with the eight global modes (the dashed cnrves), and with the whole 
spectrnm of modes (the solid cnrves). The dot-dashed cnrves represent the corresponding 
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Table 1. The eigen frequencies and overlap integrals of the global modes. 



^(k) 

^b,(k) 

Qn,{k) 

Qb,(k) 

Wi 

0.521 

0.521 

0.105 

0.1065 

W 2 

-1.11 

-1.11 

0.07 

0.072 

W 3 

-0.873 

-0.872 

0.026 

0.026 

W 4 

1.257 

1.257 

0.026 

0.026 

W 5 

1.273 

1.276 

0.024 

0.03 

We 

-1.526 

-1.525 

0.0117 

0.0127 

Wt 

-1.5 

-1.5 

0.0091 

0.0096 

Wg 

-1.353 

-1.349 

0.0066 

0.0085 


quantities calculated for the model without the phase transition. One can see that the 
curves for the model without the phase transition are very close to the curves calculated 
using only the two ’main’ global modes - the result expected from our previous discussion. 
In general, the curves calculated using the global modes give a good approximation to the 
curves calculated with help of the whole spectrum excepting the curves representing the 
transfer of energy in the inertial frame, 62 at a high rotation rate, > 10. The solid curve 
determined by the whole spectrum gives much larger energy transfer in this case. This is 
due to the influence of the modes with eigen frequencies very close to —2. It is not clear 
whether this effect is physical or it is determined by numerical inaccuracies of our method. A 
calculation with signihcantly larger number of the trial functions may resolve this question. 
However, this effect is unimportant for our purposes and will not be discussed further. 

In Figures fT^HT^ we show the distribution of Wk over the planet for the eight global modes. 
Additionally, in Figure IT^ we show the distribution of H4 for two ’non-global’ modes: Wg 
with (Tg = 0.49 and Wig with uio = —0.5. These modes have very small overlap integrals. 
In the upper panels of Figure IT^ the two ’main’ modes Wi and W 2 are shown, and in the 
lower panels of the same Figure we show the same modes calculated for the model without 
the phase transition. One can see that these modes have essentially the same structure for 
the both models. Moreover, these distributions are very similar to the distributions of the 
corresponding modes for the case of a n = 1.5 polytrope (IP). Therefore, it is natural to 
suggest that the dependence of the structure of the ’main’ modes on the structure of a 
planet is rather insignihcant (see also Figures ^1 and 1221 below). Other global modes have 
structures with amplitudes that are more concentrated towards the boundary of the planet 
when compared to the ’main’ modes. It is instructive to compare the global and non-global 
modes. One can see from Figure that the non-global modes are characterised by a large 
number of patterns with alternating signs. These patterns give compensating contributions 
to the overlap integrals thus signihcantly decreasing their values. 
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Figure 15. The overlap integrals Qn,{k) shown for the planet models with M* = IMj, and four different radii: = 1.2Rj 

(solid curve), R^ = lARj ( dotted curve), = l.SRj (dashed curve), and R* = 2Rj ( dot-dashed curve). 



Figure 16. The same as Fig. but for M* = 5Mj. The solid, dotted, dashed and dot-dashed curves are for i?* = 1.03-Rj, 
lARj, l.QRj and 2Rj, respectively. 


5.4 Eigenspectra of the planets with larger masses and radii 

In Figure HSl we show the dependency of the overlap integral Qn,{k) on the eigenfrequency 
for the planet models with planet mass equal to the mass of Jupiter with four different 
radii, i?* = lARj, l.SRj and 2Rj. The general structure of the spectrum is close to 

that of the planet with M* = IMj, and i?* = IRj, see ( Figure |HI). In all cases there are 
the two ’main’ global modes. Also there are some other ’non-standard’ global modes with 
eigenfrequencies close to the corresponding eigenfrequencies of the model with i?* = IRj. 
The ’main’ modes slightly change their position with growth of planet radius, and their 

© 2002 RAS, MNRAS 000. [TII^ 






















40 P. B. Ivanov and J. C. B. Papaloizou 



Figure 17. The dependence of ej and 62 for the planet models with Mj = IM* and i?* = 1.2Rj ( solid curves), = lARj 
( dotted curves), i?* = l.SRj ( dashed curves), and il* = 2Rj ( dot-dashed curves). The curves which always pass through 
positive values correspond to £ 2 - The curves of the same type with smaller amplitude correspond to ej and €2 calculated using 
only the two ’main’ global modes. 



Figure 18. The same as Fig. nn but for M* = 5Mj. The solid, dotted, dashed and dot-dashed curves are for R* = l.OSRj, 
ISRj and 2Rj, respectively. 


overlap integrals also slightly depend on the radius, being however rather close to each other. 
It is interesting to point out that the overlap integrals corresponding to the ’non-standard’ 
global modes are smaller than those corresponding to the model with M* = IMj, and and 
i?* = IRj. Therefore, the influence of the ’main’ global modes on the energy transfer is more 
pronounced than in that case (see Figure HTj). 
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Table 2. The binned eigen frequencies and overlap integrals of the main modes for M* = IMj. 


R* / R-pi 

1.2 

1.4 

1.8 

2 

<^6,1 

0.54 

0.58 

0.65 

0.71 

Qb^l 

0.093 

0.097 

0.117 

0.12 

^b,2 

-1.1 

-1.09 

-1.076 

-1.037 

Qb,2 

0.057 

0.06 

0.054 

0.051 


Table 3. The binned eigen frequencies and overlap integrals of the main modes for M* = 5Mj. 


1.03 1.4 1.6 2 


^b,l 

0.55 

0.53 

0.57 

0.59 

Qb,l 

0.1 

0.105 

0.098 

0.096 

<^6,2 

-1.1 

-1.1 

-1.1 

-1.1 

Qb,2 

0.061 

0.069 

0.064 

0.058 


In Figure ITT)1 we show the results of calculations for models which have M* = 5Mj with 

= 1.03-Rj, 1.6-Rj and 2Rj. As seen from this Figure the spectrum is quite stable 

for this case, such that the curves corresponding to different planet radii can be hardly 
distinguished from each other. Also, the spectrum is totally dominated by the two ’main’ 
global modes with ai ~ 0.55 and a 2 ~ —1.1. In fact, for these models it is possible to identify 
reliably only three non-standard global modes with a^k) ~ —0.85, 0.3 and 1.3. These modes, 
however, have rather small overlap integrals. 

In Fjgures IT7I and ITKI we show the dependence of the quantities 6 / and 62 which determine 
the energy transfer in the inertial and rotating frames, respectively, on for the set of planet 
model considered in this Section. Also, we compare these quantities with the same quantities 
calculated using only the two ’main’ global modes and the ’binned’ overlap integrals and 
eigenfrequencies discussed above. The corresponding eigenfrequencies and overlap integrals 
are shown in Tables |21 and El 

In Figure El we show the results of calculations for the planet models with M* = IMj 
and in Figure El the case of M* = 5Mj is shown. As seen from these Figures, the quantities 
determined by using only the two ’main’ modes give a good approximation for all curves 
excepting the cases of the planet models that have M* = IMj with i?* = 1.2Rj and 
M* = 5Mj with i?* = 1.03i?j. In those cases the curves giving 62 deviate signihcantly when 
is large. 

Figures Unil^ illustrate the effects of changing planet radius and mass on the form of 
W(k)i for the global modes. 

In Figure El we show the distributions for the two ’main’ global modes W\ and IF 2 in the 
case of the planets with M* = IMj and radii larger than the radius of Jupiter. Comparing 
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Figure 19. The distribution of for M* = IMj with two different radii i?* = lARj and R* = 2Rj. The two ’main’ 

global modes Wi ( upper panels) and W 2 ( lower panels) are shown with the left (right) plots corresponding to R* = lARj 
= 2i?j). The upper left (right) plot corresponds to fji = 0.58 with Qn^i = 0.097 (fJi = 0.71 with Qn,i = 0.11). The lower 
left (right) plot corresponds to <72 = —1.1 with Qn ,2 = 0.038 (cr 2 = —1-06 with Qn ,2 = 0.051). 

this Figure and Figure IT^ we see that the spatial structure of the ’main’ global modes only 
slightly changes with change of planet radius. These modes may be considered as ’stable’ 
modes with respect to changes of this parameter. In Figure 1201 we show the influence of 
changing radius on the ’non standard’ global mode W 3 with ~ —0.9. The distribution of 
this mode is practically the same for planet radii in the range 1.2i?j — 2Rj, and therefore 
this mode is particularly ’stable’. Examples of ’non-stable’ modes are shown in Figures IQD 
and 1221 We plot the mode IF 4 with <74 ~ 1.3 in Figure 1211 and the mode IF 5 with <75 ~ —1.5 
in Figure 1221 One can see from these Figures that the region with non-negligible hF(fc) is gets 
smaller with radius. Accordingly, the overlap integrals get smaller as well, see Figure ITHl 
As follows from Figure ITI)1 the effective spectra of the models with = 5Mj is pretty 
stable with respect to a change of radius, and the overlap integrals corresponding to the ’non 
standard’ global modes are small compared with those corresponding to the ’main’ global 
modes. Therefore, in Figure I2HI we show only the two ’main’ global modes for i?* = 1.03i?j 
and i?* = 2Rj. The general structure of the modes is similar to the case M* = IMj (see 
Figures ^1 and However, there is an interesting feature of the distributions in the 
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Figure 20. The form of the ’non standard’ global mode Wz with 0-3 ps —0.9 for the planet models with M* = IMj and four 
different radii: i?* = upper left plot; lARj^ upper right plot;i?* = lower left plot; i?* = 2 i?j, lower right plot. 


case i?* = 1.03-Rj. In this case a small scale pattern is superimposed on the large scale 
distributions of Wi and W 2 . Most probably, this effect can be explained by mixing between 
these two global modes and neighbouring local modes as a result of a perturbation. 


6 ROTATION RATE AND ENERGY EXCHANGE IN THE STATE OF 
PSEUDO-SYNCHRONISATION 

When the state of pseudo synchronisation is maintained with the value of the angular velocity 
being equal to its equilibrium value the angular momentum transfer is absent A(f2eg) = 0. 
In this case the energy transfer in the inertial frame is equal to the energy transfer in 
the rotating frame, and therefore the corresponding curves in Figures IHl El and IT^ 
intersect each other. One can see from these Figures that the intersection point have similar 
coordinates for all models we consider, and therefore the dimensionless rotation rate Ogg 
and the quantity e* determining the energy exchange through equation (m only slightly 
depend on the parameters of the planet models. We summarise results of calculations of 
these quantities in Tables El and 0 
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Figure 22. Same as Fig. 1201 but for the global mode with 0-5 ^ —1.5. 
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Figure 23. The same as Fig. I19l but for M* = 6Mj. The two ’main’ global modes Wi (upper panels) and W 2 ( lower panels) 
are shown with the left (right) plots corresponding to (i?* = 2Rj). For the upper plots cri 0.57 and Qn,i ~ 0.09 

For the lower plots (72 ~ —1.1 and Qn ,2 ~ 0.06, for the both models. 


Table 4. e* and fleg for the planet models with M* = Mj. 


R*/Rj 

1 

1.2 

1.4 

1.6 

1.8 

2 

e* 

5.3 • 10-3 
1.46 

3.6 ■ 10-3 
1.52 

3.45 ■ 10-3 
1.56 

3.6 ■ 10-3 
1.59 

3.7 ■ 10-3 
1.6 

3.6 ■ 10-3 
1.61 


As follows from these table apart from the case of a planet with Jupiter mass and radius 
which gives a rather large value e* = 5.3 ■ 10“^, the values of e* and Cleg are rather close to 
each other. Therefore, in our estimates of the time scale of tidal circularisation below we 
will use the ’typical’ values e* = 3.6 ■ 10“^ and Cleg = 1.55. Accordingly, the energy transfer 
will be given by 
3 . 6 - 10-3 


E = 

^ps 


-E^ 




(98) 


regardless of the planet radius and mass. 


Table 5. e* and for the planet models with M* = bMj. 


R*/Rj 

1.03 

1.2 

1.4 

1.6 

1.8 

2 

e* 

^eq 

3.6 • 10-3 
1.55 

3.6 ■ 10-3 
1.55 

4.2 • 10-3 
1.54 

3.7- 10-3 
1.54 

3.55 ■ 10-3 
1.54 

3.4 ■ 10-3 
1.55 
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In FigureHniwe show the dependence of Epg on r/ as a solid curve. The dotted, dashed and 
dot-dashed curves correspond to the energy transfer associated with the fundamental modes, 
Ef. They are calculated with help of equations derived in IP and presented in Appendix 
In a similar manner to the energy transfer associated with the inertial modes, Ef depends 
on the angular velocity of the planet, and we present several different curves calculated for 
different values of hi. The dashed curves are calculated for a non-rotating planet with = 0, 
they correspond to a maximal energy transfer as a function of hi. The dotted curves 

are calculated for the case hi = fleq = 1.55. The dot-dashed curves are calculated for the 
equilibrium angular velocity of the planet, associated with the fundamental modes, see 
equation dSll) of Appendix. These curves give a minimal energy transfer as a function 
of n, for a given rj. 

As follows from Appendix, in order to calculate Ef the eigenfrequencies associated with 
the fundamental modes, a/ as well as the corresponding dimensionless overlap integrals, Qf 
must be known. These quantities depend on the planet structure, and therefore they are 
different, in principle, for planets with different masses and radii, see IP for their numerical 
values. In Figure ITHl we show the case of a planet with M* = Mj and i?* = Rj (the upper 
dashed, dotted and dot-dashed curves), where we have af ~ 1.214^. and Q/ ~ 0.56. The lower 
dashed, dotted and dot-dashed curves represent the case of a planet with the same value of 
mass and i?* = 2Rj with cx/ ~ 1.712* and Qf ~ 0.44. As follows from the results provided 
in IP, when M* = 5Mj the dependence of the eigenfrequencies and overlap integrals on the 
planet radius is rather weak, and the corresponding energy transfer is always approximately 
represented by the upper dashed and dotted curves in Figure ITHl 

As follows from Figure ITCIthe inertial modes dominate the energy transfer for sufficiently 
large values of rj. When 12 = Cleg the contribution determined by the inertial waves dominates 
over that corresponding to the fundamental modes for r] > 9 in the case of i?* = Rj and 
for T] > 5.5 in the case of i?* = 2Rj. Thus, we expect that the circularisation of the planet 
orbits with sufficiently large periastron distances is mainly determined by the contribution 
associated with the inertial modes. 


® Note a misprint made in (IP). Their equations (60), (61) and (64) must by multiplied by factor tt^. This has no influence on 
the conclusions of IP. The correct expressions are shown in Appendix. 
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Figure 24. The energy transfer associated with the inertial modes (the solid curve) as well as the energy transfer associated 
with the fundamental modes as a function of the parameter 77 . The dashed, dotted and dot-dashed curves correspond to 
SI = Qeq, 0 = 0 and O = respectively. Upper (lower ) curves of the same type corresponds to a planet with M+ = Mj 
and m = Rj {R* = 2Rj). 


7 TIME SCALE OF ORBITAL CIRCULARISATION 

In this Section we apply the results developed in the previous Sections to estimate the effect 
of tides exerted by the central star on the planet on the planet orbital evolution. We assume 
that the state of pseudo-synchronisation with O = 1.55 is maintained in a situation 

when the inertial modes dominate over the fundamental modes in tidal response. 

Apart from the fact that we mainly consider the effects of inertial waves, the formula¬ 
tion of the problem is analogous to what was discussed in IP. Let us consider a planetary 
system consisting of a few massive planet and a central star and assume that a certain 
planet has lost a signihcant amount of its orbital angular momentum, most probably due to 
gravitational interaction with other planets in the planetary system. Such situations have 
been frequently observed in numerical simulations of gravitationally interacting planets (Pa- 
paloizou & Terquem 2001, Adams & Laughlin 2003). The orbital semi-major axis of this 
planet, in this case, could be of the order of its ’typical’ value 1 — lOOAu while the perias- 
tron distance Vp could be very small, of the order of a few stellar radii. For a such eccentric 
orbit the orbital angular momentum per unit of mass is related to the periastron distance 
as Lorb ~ ^2GMrp. 

If we assume that the orbital evolution is determined by tides, the angular momentum 
is approximately conserved while the semi-major axis may decrease due to transfer of the 
orbital energy into oscillations of the planet by tides and subsequent dissipation. Then, the 
energy may be radiated away from the planet. The conservation of the angular momentum 
leads to the well known fact that the radius a ’hnal’ circular orbit, rj is twice as large as 
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the periastron distance Vp. Therefore, Vp can be related to the observed orbital period of the 
planet after circularisation, Pobs = as 

rp = ^ R., (99) 

where P 3 = Pots/i^days), Rg is the stellar radius and Rq = 7 ■ 10^°cm is the radius of the 
Sun. We have, accordingly, from equation (P|) 


n = 9.ii)o ( 


Rg 


P. 


Vo = 


^ M Rr 


( 100 ) 


For the planet and the star with Jupiter and solar values of mass and radius, respectively. 


Vo ~ 1- 


Let us consider the evolution of the semi-major axis of the eccentric planet orbit due 
to tidal interactions near the orbit periastron. We can assume that the contributions from 
each of successive impulsive energy transfer can be added to each other only when one 


of two conditions is fulhlled: either the energy is dissipated during two successive perias¬ 
tron passages or the condition for the stochastic instability (|H1 is satisfied. In the latter 
case the evolution of the semi-major axis is stochastic, and we assume below that, in this 
case, our evolution equations are written for averaged over some distributions values of the 


corresponding quantities. 

Under these assumption we can write the evolution equation for semi-major axis, a{t), 
in the form (IP) 
hio _ 1 

flio 

where oio = a/(10An). The characteristic evolution time 


( 101 ) 


fio = 7.4 ■ 10^ 


' MMj 
MqM, 


R* 


10 


-9' 


^DT 


yr. 


( 102 ) 


depends on the energy transfer due to dynamic tides associated with the inertial and funda¬ 


mental modes in the planet, Epg and Ef, respectively, and also on the transfer of the orbital 
energy into stellar oscillations determined by dynamic tides exerted on the star 
{Epg + Ef + Eg) 



where Epg is given by equation and Eg stands for the energy transfer associated with 
dynamic tides in the star. Note that the time tio implicitly depends on time t through 
dependence of the planet radius P* on time due to the evolutionary cooling and a possible 
tidal heating of the planet. 

Equation (uni can be integrated to give 
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orb 


Figure 25. The circularisation time scale tev (in units of lO^yr) for IMj as a function of the planet orbital period after the 
stage of tidal circularisation, Porb (in units of days). The solid curves correspond to tev determined solely by the inertial modes. 
The dashed curves are calculated with help of the sum of contributions corresponding to the inertial modes, the fundamental 
modes and the modes excited in the star. The energy transfers determined by the inertial and fundamental modes are given 
by expressions corresponding to rotation of the planet at the equilibrium angular velocities Qps and lips, respectively. The 
dotted curves correspond to the case of a non-rotating planet. Accordingly, the contribution of inertial modes is not taken into 
account. The lower (upper) curves of the same type correspond to ain = 1 [a = lOAu] (ain = 10 [a = lOOAti]). 


®lo(0 ~ 


dt 




(104) 


where ain{tin) is the ’initial’ value of the semi-major axis expressed in units of lt]Au. A more 
accurate compared to ()102j) dehnition of the evolution time scale, tev, can be obtained from 
equation (uni requiring that the semi-major axis Oio is formally equal to zero when t = tev 
Thus 

1 /■*«’' dt 




tio{t , 


= 1 . 


(105) 


For a constant value of fio, tev ~ Ici this case we can obtain a simple estimate 

of tev assuming that the transfer of energy is determined only by the inertial modes. Setting 
Eg = Ef = 0 in equation (uni and substituting equations iHD, (unni) and (uni into (uni, 
we have 


te 


2 ■ 10 ^ 


Mj 


' M 

Mo, 


-2 


bTj 


PS 


yr. 


(106) 


Note a very strong dependence of tev on i?* and Porb- 

Results of numerical calculations of tev are shown in Figure ESI for the case of a planet 
with M* = IMj, and in Figure ESI for a planet with = 5Mj. In order to calculate tev 
we use expressions for the energy transfer associated with the fundamental modes given 
in Appendix, the corresponding eigenfrequencies and overlap integrals are given in IP. The 
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Figure 26. The same as Fip-nre l^ hiit for M* = 5Mj. 

dependence of i?* on time is taken from Burrows et al (1997) where a cooling history of 
isolated planets is considered. Thus we neglect effects of tidal heating on the planet size 
and structure reserving their discussion for a future work. Burrows et al (1997) consider 
evolution of the planet radius starting from the initial time ~ lO^yr and until a maximal 
time tmax ~ lO^yr is reached. We set the initial time Un = 2 ■ lO^yr. We assume that the 
planet radius is not changed significantly for t > tmax and set R^{t > tmax) = R*(tmax) with 
R*{tmax) = 7.2 ■ lO^cm for M* = Mj and R*{tmax) = 7.67 ■ lO^cm for M* = 5Mj. The lower 
curves of the same type in Figures EHl and EEl correspond to = IQAu and the upper curves 
are calculated for ain = lOOAu. 

In Figure ESI the solid curves describe the dependence of t^v on Porb calculated under 
assumption that only inertial waves contribute to the tidal evolution, and accordingly, in 
this case Egt = F'j = 0 in equation (uni. For the dashed curves we take into account 
all the contributions to the sum in (UnsD In order to calculate Eg we model the star as 
a non rotating n = 3 polytrope and take into account the contributions of f and low 
order p modes, see IP for details. The contribution of the fundamental modes operating in 
the planet, Ef is calculated for the case of rotation of the planet at the equilibrium angular 
velocity associated with the fundamental modes, see equation (Hi of Appendix. In this case 
Ef is given by equation (Ell). This may be justified according to the following arguments. 
The inertial modes dominates over the fundamental modes at sufficiently large values of rj, 

^ Note that conditions for the dissipation of the mode energy and the condition for the stochastic instability to set in are 
different for the inertial and fundamental modes in the planet and also for the modes determining the tidal response of the star. 
Therefore, simple adding of different contributions to the sum in is, strictly speaking, not correct. However, we assume 

hereafter that this procedure gives correct order of magnitude estimates reported in this Paper. 
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and accordingly at sufficiently large values of Port, and the fundamental modes dominate 
at smaller values of Porb- Taking into account that the contribution of fundamental modes 
decreases exponentially with increase of Porb, the transition from the case of domination of 
the fundamental modes to the case of domination of inertial modes occurs near a certain 
value of Porb- Taking into account that a planet with a low value of the moment of inertia 
has the angular velocity always close to the equilibrium value determined by the modes 
dominating the energy transfer, the dashed curve could give an order of magnitude estimate 
of tev for the whole range of Porb- The dotted curves show the results of calculation of tev for 
the case of a non rotating planet. In this case the energy transfer determined by the inertial 
modes is absent and we set Epg = 0 in equation (unsi) and we use equation ()A5fl to calculate 

Ef- 

As seen from Figure EHl when = IMj the inertial modes dominate for Porb > 4 — 
A.bdays. The time scale of circularisation is smaller than a typical life time of a planetary 
system ~ 5Gyr when Porb < Qdays for Oin = lOAu and when Porb < 5days for = 
lOOAu. It is interesting to note that the contribution of the fundamental modes sharply 
decrease when t^v > lO^yr. This is simply because the energy transfer associated with the 
fundamental modes is highly sensitive to the radius of the planet (see Appendix). It can 
lead to a substantial change of the semi-major axis only when the planet radius is close 
to its ’initial value’ ~ 2Rj. The cooling of the planet results in a decrease of the 

planet radius on a time scale ~ Itf’yr and the contribution of the fundamental modes 
quickly become unimportant after this time is reached. The tides exerted on the star are not 
signihcant at all for a planet with M* = IMj. 

In Figure |2ni we show the results corresponding to a planet with M^, = 5Mj. Contrary 
to the previous case now the inertial modes are not important for the whole range of Porb 
and the circularisation time scale at large values of Porb is determined by tides exerted 
by the planet on the star. This may be explained as follows. At large values of Porb the 
circularisation time scale is of the order of several Gyr. A planet with its mass in the range 
IMj — 5Mj and its age of the order of a few Gyr has the radius close to R*{tmax) ~ 1-Rj- The 
circularisation time scale is mainly determined by the expressions for the energy transfer 
calculated for this radius. Assuming that the inertial modes dominate we can use equation 
ffTO to estimate te^. As follows from this equation tev oc M^. When eg. Porb = Qdays and 
Oin = IQAu we have tev ~ QGyr for M* = IMj and tev ~ 25Gyr for M* = 5Mj. Obviously, 
the inertial modes cannot lead to a signihcant decrease of the semi-major distance in the 
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latter case over the life time of the planetary systems. Unlike the sitnation with the inertial 
modes, the energy transfer determined by the stellar modes grows with the planet mass, and 
therefore the tides exerted on the star dominate when the planet mass is snfficiently large. 


8 CONCLUSIONS AND DISCUSSION 

In this Paper we have developed and extended a new self-adjoint formalism for the problem 
of small stellar oscillations proposed in onr recent paper (PI). We went on to apply this to 
the calcnlation of the tidal response of nniformly rotating fnlly convective planets and stars 
which nndergo a single enconnter in a parabolic orbit or a seqnence of mnltiple enconnters 
in a weakly bonnd orbit. This was then applied to the problem of the tidal circnlarisation 
of the orbits of the extrasolar planets. 

In section 2 we began by showing that the general formalism, presented in PI when 
W = p'(?Ip + is nsed as variable characterising linear eigenmodes, can also be applied 
when the Lagrangian displacement is nsed. The advantage of this approach is that, as no low 
freqnency approximation is needed, it allows ns to consider pnlsations of arbitrary freqnency. 
We derived general expressions for the energy and angnlar momentnm transfer that occnrs 
when a rotating planet or a star passes throngh periastron in a parabolic or highly eccentric 
orbit aronnd a central mass The amplitndes and phases of the eigenmodes excited in the 
planet or a star as a resnlt of the tidal enconnter were also determined. 

We showed that these expressions can be represented in terms of a spectral decomposition 
over the normal modes of the self-adjoint operator m- We also obtained a self-consistent 
expression for the rate of dissipation of an excited mode to leading order in the small 
parameter defined by the ratio of pnlsation period to a characteristic dissipation time scale. 

In section 3 we focnsed on the case of low freqnency inertial modes, showing how the 
more general formalism in terms of the Lagrangian displacement can be rednced to that in 
PI in the low freqnency limit. This procednre allows one to separate ont the contribntion 
associated with the inertial modes and simplify the eqnations governing free pnlsation of a 
rotating planet or a star. In that case the corresponding self-adjoint problem is formnlated 
throngh eqnations and (HUll . Note that the tidal response of the inertial modes to 
impnlsive tidal forcing has not been described nsing other methods. 

We showed in section 4 that when the rotating planet or a star approaches periastron 
in an nnpertnrbed state the tidal energy and angnlar momentnm transfered throngh the 

© 2002 RAS, MNRAS Onn.lTII^ 


Dynamic tides in rotating objects 53 
excitation of inertial modes are given by very simple expressions (IHHj) and dZlD. We considered 
the mnlti-passage problem assnming that there was no signihcant dissipation of mode energy 
between snccessive periastron passages, and found a simple condition (IHH) for the occurrence 
of stochastic instability that results in the stochastic gain of inertial mode energy over many 
periastron passages. This was quantitatively similar to that given by IP for / modes. We 
found that for stochastic instability the circularisation process has to start with a >~ 30AU, 
for hnal periods of ~ 3 days, or ~ 1 — 2^417 for final orbital periods ~ 1.2 days. 

In section 5 we applied our formalism to fully convective rotating giant planets. In order 
to evaluate the expressions for the energy and angular momentum transfer the eigenspectrum 
of ()4f)|l must be found numerically, for a given model. We calculated the eigenspectrum of 
the inertial modes for several planet models with a realistic equation of state. The details 
and stability of the numerical method we used are discussed in section 5.2. We considered 
planet models with two different masses equal to IMj and 5Mj, with radii in the range 
1 — 2Rj. We found that, as in the case of n = 1.5 polytrope considered by PI, the tidal 
response was determined by a few ’global’ eigenmodes with a large scale distribution of 
perturbed quantities. Two ’main’ or ’standard’ global modes have eigenfrequencies close 
to those corresponding to the n = 1.5 polytrope. These modes have the largest overlap 
integrals characterising coupling between the tidal held and the eigenmodes. However, we 
found that for realistic planet models, there could be ’non-standard’ global modes possibly 
related to a sharp change of structure occurring near to the point of ionisation of hydrogen. 
The structure of the ’non-standard’ modes as well as their stability with respect to a change 
of the planet mass and radius was discussed. 

The results obtained for the eigen spectra of realistic planet models were applied to the 
problem of tidal circularisation of the orbits of the extrasolar planets in sections 6 and 7. In 
section 6 we calculated the rotational angular velocity for which the net angular momentum 
transferred was zero. Because of the relative low inertia of the planet compared to the 
orbit it is expected to rapidly attain this angular velocity and achieve a state of pseudo 
synchronisation. The angular velocity associated with pseudo synchronisation was found to 
be always close to 1.55 times the angular velocity of a circular orbit at periastron. 

Orbital evolution arises through the transfer of orbital energy to the modes of oscillation. 
Transfers at successive periastron passages lead to orbital circularisation. This is the case, 
either when the mode energy is dissipated directly between encounters, or when the system 
is in the regime of stochastic instability. We compared the contribution associated with the 
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inertial modes with the contribution associated with the fundamental modes, in a state of 
pseudo synchronisation, to the circularisation time scale in section 7. We found that the 
inertial modes led to effective circularisation of a giant planet with mass ~ IMj on a time 
scale smaller than or of the order of a few Gyr, when the initial semi major axis was less 
than 1{)AU and the final orbital period after circularisation Porb < Qdays. In this case the 
inertial modes play the most important role at the longer periods. However, the inertial 
waves are less important for the planets of larger mass 5Mj that we considered. In that case 
the dynamic tides exerted on the star central star, which we assumed to be solar like, play 
a major role. 

In our opinion, the most important unresolved question related to the theory discussed 
above is concerned with internal dissipation. Of particular concern is the internal location of 
modal dissipation and the related issue of whether the energy deposited can be radiated away 
between successive periastron passages. There are several potentially interesting channels of 
dissipation which may operate in our case. The non-linear parametric instability discussed 
by Kumar and Goodman (1996) for the case of g modes in convectively stable stars may be 
effective for the inertial modes as well. This instability operates when there are eigenmodes in 
the stars with eigenfrequencies approximately half that of the eigenfrequencies of the tidally 
excited modes and sufficiently large coupling coefficients. These modes may be unstable 
on account of parametric resonance. Taking into account that the eigenspectrum of the 
inertial modes is dense, unstable modes are expected to exist for sufficiently weak dissipative 
processes. This question will be addressed in our future work. 

Our formalism is general enough that it can be extended to the case of convectively stable 
uniformly rotating stars and it can also be applied to other astrophysical situations where 
excitation of stellar pulsations with frequencies comparable to the angular velocity of rotation 
is important. In this particular work we assumed that the orbital and planetary angular 
momentum vectors were aligned. This is a very reasonable assumption when the system has 
many periastron passages and so readily attains a state of pseudo synchronisation during 
the circularisation process. However, misaligned angular momenta should be considered if 
the main issue is tidal capture, as in globular clusters, because in that case there is no reason 
for the initial angular momentum vectors to be correlated. These issues and extensions of 
the formalism will be considered in future work. 
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APPENDIX A: ENERGY AND ANGULAR MOMENTUM TRANSFER 
ASSOCIATED WITH THE FUNDAMENTAL MODES 

Here we show the expressions for energy and angular momentum transfer associated with 
the the fundamental modes in the limit r/ —> cx) (IP). For the energy transfer we have 



(Al) 


and for the angular momentum transfer we have 



(A2) 


where 



(A3) 
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and 


F2{n) = 1 


9 


16^2 
re 3 


Pflr] 


(A4) 


2^^{afr])^ 

df = af/Q^ and f2 = are the dimensionless eigenfreqnency and dimensionless angular 

velocity of the planet, respectively and a)+ = d/ + 2/3f2. The dimensionless quantity (3 
determines the frequency splitting of the mode due to rotation, it is close to 0.5 for our 
models (see IP and references therein). Qf are the dimensionless overlap integrals, see PT 
and IP. When = 0 the energy transfer as a function of is maximal. 


rpmax 

A/ 


16^2 




-TT a 


fQM 


(A5) 


Assuming that the fundamental modes dominate the tidal response, the state of pseudo¬ 
synchronisation is determined by condition = 0. In this case we obtain from equation 

dH 

nL = 


ps 


In I 




4V2f3r] 

When the energy transfer is minimal 


Tpmin 

A/ 


5v/2 T] 




(A6) 


(A7) 
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